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Abstract. We investigate the influence of aperiodic perturbations on the critical 
behaviour at a second order phase transition. The bond and site problems 
are compared for layered systems and aperiodic sequences generated through 
substitution. In the bond problem, the interactions between the layers are distributed 
according to an aperiodic sequence whereas in the site problem, the layers themselves 
follow the sequence. A relevance-irrelevance criterion introduced by Luck for the bond 
problem is extended to discuss the site problem. It involves a wandering exponent 
for pairs, which can be larger than the one considered before in the bond problem. 
The surface magnetization of the layered two-dimensional Ising model is obtained, 
in the extreme anisotropic limit, for the period— doubling and Thue-Morse sequences. 



1. Introduction 

The discovery of quasicrystals (Shechtman et al 1984) has opened a new field of 
research which has been quite active during the last ten years (see Henley 1897, Janssen 
1988, Janot et al 1989, Guyot et al 1991, Steinhardt and DiVicenzo 1991). On the 
theoretical side quasiperiodic or, more generally, aperiodic systems are interesting 
because they appear as intermediates between periodic and random ones. Thus, 
phase transitions in such systems are expected to display a rich and unusual critical 
behaviour. 

Studies of the Ising model (Godreche et al 1986, Okabe and Niizeki 1988, S0rensen 
et al 1991), the percolation problem (Sakamoto et al 1989, Zhang and De'Bell 1993) 
and the statistics of self-avoiding walks (Langie and Igloi 1992) on the two-dimensional 
Penrose lattice did not show any change of critical exponents. Universal behaviour was 
obtained in three dimensions too (Okabe and Niizeki 1990). But some systems were 
also found for which the aperiodicity has some influence. One may mention interface 
roughening in two dimensions for which a continuously varying roughness exponent 
was obtained with the Fibonacci sequence (Henley and Lipowsky 1987, Garg and 
Levine 1987). 

Some exact results have been obtained for the layered Ising model with an aperiodic 
modulation of the interlayer couplings (Igloi 1988, Doria and Satija 1988, Benza 1989, 
Henkel and Patkos 1992, Lin and Tao 1992, Turban and Berche 1993). The problem 
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was studied in the extreme anisotropic limit where it leads to a one-dimensional 
aperiodic quantum Ising model (QIM) in a transverse field, which is often easier 
to handle (Kogut 1979). For the Fibonacci and other sequences, the specific heat 
was found to display the Onsager logarithmic singularity but for different sequences 
the singularity is washed out (Tracy 1988), like in the randomly layered Ising model 
(McCoy and Wu 1968a, 1968b, McCoy 1970). 

The situation was clarified in a recent study of the bulk properties of the QIM 
(Luck 1993a). In this work. Luck proposed a generalization for aperiodic systems of 
the Harris criterion (Harris 1974), allowing a classification of critical aperiodic systems 
according to the sign of a crossover exponent This exponent involves the correlation 
length exponent of the unperturbed system v and the wandering exponent lj governing 
the size-dependence of the fluctuations of the aperiodic interactions (Dumont 1990). 
The criterion has been also applied to the case of anisotropic critical systems with 
uniaxial aperiodicity (Igloi 1993), explaining the interface roughening results. It was 
later generalized to d-dimensional aperiodicities in isotropic critical systems (Luck 
1993b). The form of the singularities with a relevant aperiodic perturbation has been 
discussed by Igloi using scaling arguments (Igloi 1993). Recently, some exact results for 
the surface magnetization of the QIM have also been obtained for irrelevant, marginal 
and relevant aperiodicities (Turban et al 1994, Igloi and Turban 1994) and conformal 
aspects have been discussed (Grimm and Baake 1994). 

Most of the systems treated so far were dealing with an aperiodic distribution of 
the couplings, i.e. with the bond problem. In magnetic systems, this corresponds to 
an aperiodic distribution of the atoms mediating the interactions in a supcrcxchange 
mechanism. In the present work, we study the surface magnetization of the aperiodic 
QIM, comparing the bond problem examined previously (Turban et al 1994, Igloi and 
Turban 1994) to the site problem for which the magnetic moments are distributed 
aperiodically and interact through a direct exchange mechanism. Then the couplings 
depend on the nature of the two atoms involved in the interaction. We show that, for 
a given aperiodic sequence, the perturbation may be more efficient for the site than 
for the bond problem and may lead to a different critical behaviour. Exact results are 
obtained for two typical aperiodic sequences. 

In section 2 we present the Hamiltonian of the QIM and give the expression of 
the surface magnetization, defining the parameters for the bond and site problems. 
In section 3 we recall the properties of the substitution matrix, associated with an 
aperiodic sequence generated through an inflation rule, for the bond problem. The 
substition matrix adapted to the site problem is defined and compared to the previous 
one in section 4. Then the relevance-irrelevance criterion is discussed (section 5) and 
some general results about the QIM critical coupling and surface magnetization are 
presented (section 6). The period-doubling and Thue-Morse sequences are studied in 
sections 7 and 8 and the final section contains a summary and discussion. 



2. Hamiltonian and surface magnetization 

Let us consider a layered semi-infinite two-dimensional Ising model with exchange 
interactions K-y{k) parallel to the surface and K2{k) between the layers at k and 
fc + 1 (in k^T units). The extreme anisotropic limit (Kogut 1979) corresponds to 
K-^{k) — > oo, K2{k) while keeping the ratio = K2{k)/ Kl{k) fixed. In this 
expression Kl{k) — —1/2 In tanh if ^ (A:) is a dual coupling which goes to zero in the 
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limit. Introducing a constant reference coupling K^, the dual coupling can be rewritten 

as h^,Kl where h^, is the transverse field. The transfer operator QyiY>{—2K11-C) involves 
the Hamiltonian of a one-dimensional spin 1/2 quantum chain. Introducing the two- 
spin interactions Jj. = /i^A;,, the QIM takes the following form 

oo 

^ = -2E['^^^^ + 4^fcX+i] (2.1) 
^ fc=i 

where the crs are Pauli spin operators. 

The Hamiltonian can be put in diagonal form (Lieb et al 1961) 

n = Y.eJvl%-l\ (2.2) 

using the Jordan-Wigner transformation (Jordan and Wigner 1928) followed by a 
canonical transformation to the diagonal fermion operators r]^, . The fermion excitation 
spectrum is obtained as the solution of the eigenvalue problem 

^a^aik) = -hMk) - JMk + 1) 

eMk) = -Jk-Mk - 1) - h^k) (2.3) 
/iq = Jq = 

where the (/)„(fc) and ipaik) are the components of two normalized eigenvectors which 
satisfy the boundary conditions (/)^(0) =■0^(0) =0. 

In the ordered phase, the two spin correlation function in the surface asymptoti- 
cally gives the square of the surface magnetization, which can be written as the matrix 
element =< llaJ'lO > between the ground state and the first excited state of the 
Hamiltonian. For the semi-infinite system these two states become degenerate in the 
ordered phase, i.e. the lowest excitation vanishes. Using the above-mentioned 
transformation to diagonal fermions, it can be shown that is also given by 4>i{i)- 
According to the first equation in (2.3) with = 0, the other components of the 
eigenvector can be deduced from the recursion relation 

4,^{k+l) = -^Mk) = -Xk'Mk). (2.4) 

The normalization of the eigenvector then leads to the surface magnetization 
(Peschel 1984) 

\ -1/2 
oc J \ 

i+En\" ' (2-5) 
j=i fe=i / 

where the couplings Jj, and hf, in Ti. only enter through their ratio A^,. Thus is the 
same as for a quantum chain with h/^ = 1 and an effective two-spin interaction A;;,. 
Such a reparametrization of the Hamiltonian is no longer possible when nonvanishing 
excitations are involved, if the transverse field is k dependent. It follows that, in 
general, the effect of both interactions have to be considered. We shall come back to 
this point in section 5. 

In the bond problem the interactions parallel to the surface are constant, K^(k) = 
K^, so that /ifc = 1 and A^ = K2{k)/K^ depends on the layer index k only through 
the value of the interlayer interaction. In the site problem, the effective coupling 
Aj, = K2{k) / K1{k) involves both the interaction inside layer k and the interaction 
between layers k and k + 1. As a consequence, its value depends on the nature of the 
two layers and it is generally asymmetric. 
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3. Substitution matrix for the bond problem 

In this section we give a brief summary of the properties of aperiodic sequences 
generated through an inflation rule. To simplify the presentation, we consider 
sequences involving only two letters A and B. In the bond problem these letters 
correspond to the interlayer interactions A^, A^. The generalization to any number 
of letters is straightforward. 

A sequence is constructed through iterated substitutions on the two letters, 
A — > >S{^}, B — > S{B}. The process will be illustrated on the following example: 

S{A} ^BAA, 

(3.1) 

S{B} =AB. 

When the construction starts on A, after n steps, one obtains: 

n = A 
n=l BAA 

(3.2) 

n=2 ABBAABAA 



whereas, starting on B, the iteration gives 

n = B 
n = l AB 

(3.3) 

n = 2 B AAAB 



Informations about the sequence arc contained in the substitution matrix 

/S{A} S{B}. / X 

where the matrix elements give the numbers of ^ or B in S{A} or S{B}. It is easy to 
check that the numbers and of A and B in the sequence, after n substitutions, 
are given by the corresponding matrix elements in Mj". They belong to the first 
(second) column when the construction starts on A (B). 

Let be the right eigenvectors and the eigenvalues of the substitution matrix 
such that 

MiV„ = A„V„. (3.5) 

L^, and the length of the sequence, are asymptotically proportional to 
A" where A^ > 1 is the eigenvalue of Mj with largest modulus, which is real 
and positive according to the Perron-Frobenius theorem. The asymptotic density, 
pd, = lim„^^ L^/L„, can be deduced from the associated eigenvector with 



Surface magnetization of aperiodic Ising systems 



6353 



The interactions in the bond problem can be rewritten as 

Aa = A + p^.5i, As=A-p^(5i, (3.7) 

where A is the averaged couphng and = Ag — A ,^. the amphtudc of the aperiodic 
modulation. At a length scale the cumulated deviation from the average is 

in 

Ai(iJ = E (^fe -^) - -^1^2 - W (3-8) 
fe=i 

where A2 is the second largest eigenvalue in modulus. The wandering exponent 



InlA^I 



(3.9) 



governs the behaviour of the fluctuations of the interlayer couplings (Dumont 1990). 



4. Substitution matrix for the site problem 

Let us now consider the site problem. The two letters A and B then correspond to 
the two magnetic species which are distributed according to the aperiodic sequence, 
with the layers containing either A or B atoms. The QIM effective couplings A^^, 
^AB 1 ^BA \bb depend on the nature of the two layers involved in the interaction 
through the intra- and interlayer couplings in the anisotropic classical system. 

In order to count the numbers of bonds of different types in the sequence after 
n substitutions, L-^"^, L^^, L^^ and , one defines the two-letter substitution 
matrix 



/S{A[A]} 
I ^AA 



S{AIB]} S{BIA]} S{BIB]}- 



"-AA 



''AA 



'AA 



V 



S{A[A]} S{A[B]} S{BIA]} S{BIB]} 

''AB "'AB "-AB "-AB 

S{A[A]} S{A[B]} S{B[A]} S{B[B]} 

''BA "'BA "bA "-BA 



S{A[A]} S{A[B]} S{B[A\} S{B[B\} 



''BB 



''BB 



''BB 



''BB 



/I 2 0\ 
10 11 

110 1 

Vo 1 0/ 



(4.1) 



where for example, ^i^^'^'"*^ gives the number of A-B-bonds in the sequence generated 
by S{B} complemented by the first letter in Such a sequence builds the first 

part of the sequence which is obtained when the inflation rule is applied to a BA- 
bond. The same matrix Mj has been considered before as resulting from substitutions 
on words of length two (Queffelec 1987). 

With the example of the previous section, the two-letter substitutions read 



S{A[A]} = B A A [B] 
S{A[B]} = B A A [A] 
S{B[A]} = AB [B] 
S{B[B]} = AB [A] 



(4.2) 



6354 L Turban et al 



and lead to the last matrix in equation (4.1). 

As before, the matrix elements in eaeh eolumn of M2" give the mimbers of bonds 
of each type (LJ^; i,j = A,B) in the sequence obtained after n iterations. These 
numbers are found in the two first (last) columns if the construction starts on A (B) 
and are given by the minimum of the two values appearing in each half- row. Due 
to end effects, one of the values in each column differs by 1 from the true number of 
bonds in the sequence. At n = 2, for example, we have 




= 1- (4-3) 



The first two columns correspond to a sequence constructed on A at the second 
iteration and lead to L2^ = 2, L^^ =2, ^^ = 2 and ^ = 1, in agreement with (3.2). 

Since the sum of the numbers of bonds starting with a given letter (i.e. with A, 
AA and AB bonds) gives the number of times this letter is met in the sequence for 
the bond problem, the matrix elements of Mj are recovered by taking the sum of the 
two first elements and the sum of the two last elements in each column of M2. The 
same results are obtained with the two first (last) columns since the corresponding 
sequences only differ through their last bond. The same relation exists between the 
elements of Mj" and Mj"". 

Let be the eigenvalues and the right eigenvectors of Mj. The numbers of 
bonds {i,j = A,B) in the sequence after n iterations are still proportional to the 
nth power of the largest eigenvalue . Using the associated eigenvector, one obtains 
the asymptotic bond densities 

Poo ^4 TiA f ■\ v^4 



(4.4) 

BA^ W^l(3) BB _ W,{A) 

Since the length of the sequence after n steps is also the sum of the L^, the 
leading eigenvalues of the two matrices are the same. Using the above-mentioned 
relation between the matrix elements of M2 and M^, the secular equation of Mj can 
be factorized. The first factor gives back the secular equation of Mj so that fij^ = 
and A2 also belongs to the spectrum of Mj . The two last eigenvalues of follow 
from the second factor and read 

1 



« 2 



b ± V(a - hf + 4cd (4.5) 



where 



„ _ S{A[A]} _ S{A[B]} , _ S{B[A]} _ S{B[B]} 

" — ^AA '''AA " — "-BA "'BA 

„_S{A[A\} S{A[B]} , _ S{B[A]} _ S{B[B\} 

^ — '''BA "'BA " — '''AA "aA 



(4.6) 



Since the coefficients in (4.6) involve differences between the numbers of bonds in 
sequences which at most differ through their last bond, they are equal to or ±1. They 
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are completely determined through the first and last letters in S{A} and S{B} and 
can be obtained by inspc;c;tion. The two eigenvalues are also equal to or ±1. When 
the two substitutions begin with the same letter, the coefficients and the eigenvalues 
f2„ vanish. Other cases are listed in table 1. 



Table 1. Coefficients of the secular equation (4.5) and corresponding eigenvalues for 
substitutions starting with different letters. 



S{A} 


A... A 


A.. .A 


A. 


..B 


A... 


. B 


B .. 


.A 


B. 


..A 


B . 


..B 


B ...B 


S{B} 


B...A 


B...B 


B . 


..A 


B.. 


.B 


A.. 


.A 


A. 


..B 


A.. 


..A 


A...B 


a 


1 


1 












- 1 




- 1 












b 





1 







1 









- 1 









- 1 


c 


1 





1 









- 1 









- 1 







d 








1 




1 














- 1 




- 1 




1. 




1. - 


-1 


1. 




- J: 




- 1. - 


-i 


L - 


-I 


- 1. 1) 



Like in the bond problem, the fluctuations in the couplings A^, can be deduced 
from the substitution matrix, working in the basis of the right eigenvectors W^,. The 
averaged coupling A is now given by 

A = pt^'XAA + pi^'^AB + pZ^'^ba + pZ^'^bb- (4.7) 

The cumulated deviation from A is generally governed by the second largest 
eigenvalue of Mj in modulus with, at a length scale L^, 

A2(iJ = E (^fc - ^) ^^^^ ^^^n- (4-8) 

fe=l 

Here, the amplitude of the perturbation 62 is a linear combination of coupling 
differences, A^^ — A^g, A^^ — A^^, . . ., and the wandering exponent for bonds can 
be written as 

As a consequence, the wandering exponent changes when |fi2l = 1 > l-'^2l- Some 
examples are listed in the next section. 



5. Relevance— irrelevance criterion 



In both problems the aperiodic modulation introduces a thermal perturbation above 
A which, at a length scale i, has an averaged density 

JXIL) = ^ ~ 5,L-^-' [i = 1, 2) (5.1) 

where the w^s are the wandering exponents defined in equation (3.9) for the bond 
problem and equation (4.9) for the site problem. 
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Changing the length scale by h = L/L' leads to the renormalized density 

^i{L')-8[(^^y =6V''<5.i--i (i = l,2), (5.2) 

where v is the bulk correlation length exponent. As a consequence, the perturbation 
amplitude 5^ transforms like 

5[ = (5.3) 

with a crossover exponent $ • = 1 + z/(a;^ — 1) (Luck 1993a, Igloi 1993). 

When $j is positive, the perturbation grows under rescaling, leading to a new 
fixed point with a different critical behaviour. When is negative, the perturbation 
decays under rescaling and the critical behaviour is the same as for the homogeneous 
system. The crossover exponent vanishes in the marginal case and then the critical 
exponents may vary continuously with the perturbation amplitude, i.e. the system 
may display a nonuniversal critical behaviour. 

For the site problem, as mentioned at the end of section 2, the aperiodic modulation 
does not generally combine into a single effective parameter and one has to consider its 
effect on the transverse field and the two-spin interaction separately. The transverse 
field may take the values ft,^ or like Aj, in the bond problem. Thus, the behaviour 
of the corresponding perturbation under rescaling is controlled by whereas the 
perturbation of the two-spin interaction Jj, is governed by $2- Since, according to the 
discussion of section 4, lo2'>lOi, for any singular quantity the critical behaviour will 
depend on the sign of $2 i'^ the site problem, like for the surface magnetization. 



Table 2. Comparison between the bond and site problems for some typical 
sequences. The last two lines refer to the 2d Ising model. 



sequence 


Fibonacci" 


Thue-Morse'' 


period- 
doubling° 


three- 
folding'* 


fivefold*^ 


S{A} 


B 


AB 


BB 


ABA 


AAAB 


S{B} 


BA 


BA 


BA 


ABB 


BBA 


Al 


rf 


2 


2 


3 


2 + T 


IA2I 


T - 1 





1 


1 


3 - T 




T - 1 


1 


1 


1 


3-T 




- 1 


— 00 








0.25157 




- 1 











0.25157 


bond problem 


irrelevant 


irrelevant 


marginal 


marginal 


relevant 


site problem 


irrelevant 


marginal 


marginal 


marginal 


relevant 



" See e.g. (Tracy 1988). 

^ See e.g. (Dekking et al 1983a). 

This sequence appears in connection with the period-doubling cascade (Collet and Eckmann 1980). 

This is the folding sequence of a dragon curve (Dekking et al 1983b). 

This sequence is connected with tilings of the plane with fivefold symmetry (Godreche and Luck 
1992, Godreche and Langon 1992). 
t=(1+a/5)/2 is the golden mean. 

Table 2 gives the two largest eigenvalues of the substitution matrices and the 
wandering exponents for typical aperiodic sequences. For the two-dimensional Ising 
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model with i/ = 1, the borderline between relevant and irrelevant behaviour corresponds 
to ujj ={) according to (5.3). The bond and site aperiodic perturbations arc irrelevant 
for the Fibonacci sequence whereas the site perturbation becomes marginal for the 
Thue-Morse sequence. The following sequences all have divergent fluctuations so that 
the perturbation behaves in the same way for both problems. The period doubling and 
three-folding sequences lead to marginal perturbations but the nonuniversal exponents 
should differ for the bond and site problems. The last sequence is relevant in both 
cases. 



6. General results: critical coupling and surface magnetization 

The critical coupling of the inhomogeneous QIM is such that (Pfeuty 1979) 

.'»fi(i^)"'=,'™riw:"='- 

fe=l ^ k/ c 

For the bond problem, we take = A as a refcrence coupling and A^ = Ar. We 
associate a digit with the A:th letter in the sequence. With /fc = for A and /fc = l 
for B, the A;th coupling can be written as A;. = Xr^*" . The critical coupling A^ is such 
that lim^^^ A^r"^/^ = 1 with 



E/^- (6.2) 

k=l 

Finally, one obtains 

A^ = r-P'" (6.3) 

where — defined in equation (3.6). 

For the site problem, we take the following parametrization: 

A^^ = A, Xab = ^BA = ^BB = ^'T- (6-4) 

Due to the transverse field contribution, A^^ is generally different from A^^. The 
effective coupling can be written as 

— ^^fkfk + lj^fk + l—fkfk + lyfk—fkfk + l 

and using (6.1), the critical coupling is such that 

lim u(/^+i-/i)/-f^A,s"^/^ /r\™i./i ^ ^ 
where s = uv and 

mj = i2Mk+i- (6-7) 
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The critical coupling is then given by 



nil 
oo L 



(6.8) 



Alternatively, using the asymptotic densities defined in (4.4), the critical coupling can 
be expressed as 



(6.9) 



According to equation (2.5) and making use of (6.2) with tiq = 0, the surface 
magnetization for the bond problem is given by 



m^ = [S{X,r)] 



-1/2 



5(A,r) = ^A-2V 

3=0 



■2n,- 



(6.10) 



For the site problem, using equations (2.5), (6.2), (6.5), and (6.7) with mg = 0, we 

have 



"^s = P(A,r,s,u)] , 

^ —2 

S(A, r, s, u) = u'^f' ^" 



2nj -2/j+i 



3=0 

Since /j. = 0, 1, one may use the identity 

= 1 + (a - 1)4 
to rewrite the sum as 

S(A, r, s, u) = u^^' 

where 



(6.11) 



(6.12) 



El (a, ^, s) + (u-2 - 1)S2 (a, ^, s)] , (6.13) 



Si(A,a;,y)=^A-2^a;-2™^2/-2"^- 



j=0 

OO 



S2(A,a;,y) = ^A-2^': 

j=0 



-2mT . — 2n 



2/ Jj+1- 



(6.14) 



Let us consider the case s = r, i.e. A^^A^^ = A^^A^g, a common approximation 
in the case of symmetric couplings (A^^ = A^^). Then Si(A, 1, r) = S{X, r) and, using 
the identity fjj^i = {r~'^^'*^ — l)/(r~^ — 1), one obtains 



E(A,r,r,w)=u2(/i-i)L^ 



1-H 



AV + 



1-^2 



S(A,r)-AV2 



(6.15) 



It follows that, for this particular combination of couplings, the critical behaviour is 
governed by ^(A, r) and is the same as in the bond problem. 
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7. Period— doubling sequence 



The results of the preceding sections will be now illustrated on the examples of 
two sequences with different surface magnetization exponents for the bond and site 
problems. Wc begin with the period doubling sequence (Luck 1993a, Collet and 
Eckmann 1980) which is marginal for both problems, according to table 2. 

Since the solution of the bond problem can be found elsewhere (Turban et al 
1994), wc only give a summary of the results. Starting the iteration on B and using 
the substitutions given in table 2, one obtains 

B 

B A 

BABE 

(7.1) 

BABBBABA 

BABBBABABABBBABB 



The asymptotic density Poo = 2/3 leads to the critical coupling 

= r~^/^ (bond problem). 
The form of the substitution is such that 



/2fe — 1 //s' 

from which one deduces 



2/s+l 



2.7 — n,, n 



2j+l 



2j + l 



J' 



(7.2) 
(7.3) 

(7.4) 



Splitting the sum S{X,r) in equation (6.10) into even and odd parts and using (7.4), 
one obtains the recursion relation 



5(A,r) 



1 + 



1 



5(AV,r-i) 



and the series can be written as an infinite product (Turban et al 1994) 



5(A,r) = n 



fc=i 



1 + A, 



1 + Ar^ 



(7.5) 



(7.6) 



Let S{z) be the series expansion of S{X,r) in powers of z = (A^/A)^. According to 
equation (6.10), near the critical point S{z) should display a power law singularity 
with S{z) ~ (1 — z)~'^^% where (3^ is the surface magnetization exponent. It may be 
shown that, at the critical point z = l, the truncated series S^{z) containing the first 
L terms in S{z) behaves as L"^^' (Igloi 1986). Since the first I terms in (7.6) just 
contain the truncated series with L = 2^', the surface magnetization exponent is given 
by (Turban et al 1994) 



^ ln[(l + AJ(l + A-i)] ^ 1 ln(AyVA£^ 
'^^ 41n2 2 ln2 



(7.7) 
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Figure 1 . Temperature dependence of the surface magnetization (period-doubling, 
site problem) for different values of r, u=.b and v = 2. 
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Figure 2. Temperature dependence of the surface magnetization (period-doubling, 
site problem) for different values of u, r = .b and s = l. 



Let US now consider the site problem. The critical coupling, 
= {rs)~^^^ (site problem), 



(7.8) 



follows from (4.4) and (6.9) with p^g= p^^— pgg = l/3. Putting (7.3) into (6.7), one 
obtains the recursion relations 



m 



2i 



2j - 2n^., =2j + l- 2nj - 



(7.9) 



which can be used, together with those given in (7.4), to relate S^, T,^ and S. Splitting, 
as above, the sums in equation (6.14) into odd and even values of j, and using the 
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Figure 3. Variation of the surface magnetization exponent with r^/s in the site 
problem for the period-doubling sequence. 



identity 



^fk 



b — a a — 1 
6-1 ^ b-1 



(7.10) 



leads to: 



E,(A, X, y)=a^ + k,S [{Xxyf, {x^y)-'] , i=l,2, 
= -{Xxy) 



a;4y2 _ 1 

■2 



x^y^ — 1 



{Xxyf 

A-2x2 - {Xxyf 



(7.11) 



1 + 



a;4y2 _ 1 



These relations, with x = r/s, y = s, and equation (6.13), with = 1, finally give 



E(A, r,s,u) = a + bS 



\2 2 
AT 



\2 2 

a = A r 



2^,2 



6=1 + 



s — r u 



(7.12) 



X-\r' - u') + X'r^r'u' - s') 



which, together with (7.6), solves the site problem. Some examples of the temperature 
variation of the surface magnetization are shown in figures 1 and 2. Similar curves, 
for the bond problem, can be found in Turban et al (Turban et al 1994). 
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The critical behaviour of the surface magnetization for the site problem is governed 
by the singularity of S(\' ,r') in (7.12), at given by (7.2), with A' = A^r^, 
r' = sr~^. It follows that the critical coupling satisfies A^r^ = (sr~^)~^/^, in agreement 
with (7.8). Changing A^ into A^ = (sr~^)~^/^ in equation (7.7), one obtains the surface 
magnetization exponent for the site problem 



lln 



^^=2^ h^2 (^-1^) 

When r = s, the critical exponent is the same as for the bond problem in (7.7), as 
shown in section 6. The variation with r^/s is shown in figure 3. One may notice 
that the value = 1/2, correspondind to the surface exponent of the homogeneous 
two-dimensional Ising model at the ordinary surface transition, is recovered when 
s = uv = and A^, = r^^. This result is linked to the absence of AA pairs in 
the sequence. This value is also the minimum value of the exponent (3^, i.e. the 
aperiodicity generally weakens the singularity. 



8. Thue— Morse sequence 

As a second example, we consider the binary Thue-Morse sequence (Dekking et al 
1983a), which, as mentioned in table 2, leads to an irrelevant perturbation for the 
bond problem, treated in Turban et al (Turban et al 1994), and to a marginal one for 
the site problem. Starting with the letter A, the Thue-Morse subsitution in table 2 
give, successively, 

A 

A B 

ABBA 

(8.1) 

ABBABAAB 

ABBABAABBAABABBA 



In the bond problem, the critical coupling 

A^ = (]3on(j problem) (8.2) 

follows from the asymptotic density = 1/2 given by equation (3.6). The form of 
the substitutions immediately leads to the recursion relations 

/2fe = 1 ~ /fe' /2fe+i — /fe+i' i^-'^) 
which can be used in (6.2) to give 

= = j + fj+i- (8-4) 

A chain with length L = 2j has a density equal to the asymptotic one, which explains 
the vanishing second eigenvalue A2 in table 2. 
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Since the calculation of has aheady been described elsewhere (Turban et al 
1994), we just mention the results. The surface magnetization follows from 



oo 



(8.5) 



fc=0 



where Srpjy^(x) = Y^'^=i fk^ is the Thue-Morse series (Dekking et al 1983a). Near the 
critical point, one obtains 



1 /1-A? 



4 V 1 + A2 



t + 0{t^ 



(8.6) 



where t = 1— (A^/ A)^ is the deviation from the critical point. The surface magnetization 
exponent takes its unperturbed value, /3g = 1/2, as expected for an irrelevant 
perturbation. 

In the site problem, with = = 1/6 and = p^ = 1/3, the critical 
coupling takes the form 



A^ = r-^/'^s"^/^ (site problem). 
Making use of (8.3), equation (6.7) leads to 



(8.7) 



(8.8) 



The same result is obtained for even and odd terms because the two letters occur 
within successive pairs along the sequence, so that /2fc+i/2fe+2 = 0- 

With a sequence starting on A, the front factor disappears in the sum E(A, r, s, u) 
of equation (6.13). The S^s, defined in (6.14), satisfy the functional equations 

S.(A,a;,t/) = ai^i:.{\^y,x-^,x), i = 1,2, 



Hi 



-2 \ -2 



12 = a; - A + {XxyY 

= A-^ 



(8.9) 



21 



-^22 



X 



which, as usual, are obtained by splitting the sums into even and odd parts and 
using (8.4), (8.8), as well as the identity (6.12). 



Through iteration, at step fc > 1, the arguments of the S^s become 

3 2'=-^ -1 ( 

y ' Xk = yk =^ 



X = ;^2*'^[2*^-i + (-l)*']/3 2"-! ^ ^ , -1 ^ ^(-1)'= 



5.10) 
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Figure 4. Temperature dependence of the surface magnetization (Thue— Morse, site 
problem) for different values of r, u=.5 and s = l. 




Figure 5. Temperature dependence of the surface magnetization (Thue— Morse, site 
problem) for different values of u, r = .5 and s = l. 



and the functional equations (8.9) can be rewritten in matrix form as 

^2('^fe>^fe!2/fc)/ V^2('^/c+i>2;/j+i,yfe+i)/ ' 



with 



Then, 



1 + a: 



fc > 1. 



Si(^i.a;i,yi) 
S2(Ai,a;i,yi) 



fe=i 



(8.11) 



(8.12) 



(8.13) 
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Figure 6. Variation of the surface magnetization exponent with r/s in the site 
problem for the Thue-Morse sequence. 



where the components of the vector on the right follows from the form of the E^s when 

Equations (6.11), (6.13), and (8.11-13), formally solve the problem. The tempe- 
rature variation of the surface magnetization is shown in figures 4 and 5. 

It seems difficult to obtain an explicit expression for the surface magnetization 
since the form of the matrix Ty^ depends on the index k. But some progress can 
be made concerning the critical behaviour, using the same scaling method as for the 
period doubling sequence, i.e., looking at the L-dependence of the truncated series 
'^l{z) at the critical point z = {X^/Xf = 1. 

With the actual values of the arguments in (6.13), x = r/s and y = s, 
equation (8.10) gives 



r\(-i)V3 



Vk 



(-1)' 



l.U) 



so that, at the critical point, the matrix elements in T|^ depend on k only through 
{—!)''■ Let us introduce the product U = T2p_i Tjp with 



U 



2 + w + (w ^-w){l+w + '. 
w + w'^ (l-w)(2 + w) 



w 



2/3 



1.15) 



Taking i = 2^', the first L terms in T,j^{z) are obtained, up to L-independent factors, 

keeping the first 21 terms in the infinite product of equation (8.13). Due to the form of 
the matrix elements in (8.12), they belong to the first colmim of the matrix resulting 

from the finite product. Thus, using (8.15), they are also given by U' 

critical point. After diagonalization, one obtains 



S^(l)~(22')2A~(«,V4 + ^-l/4y 



(8.16) 
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where the last term is the Ith. power of the largest eigenvalue of U. Finally, the surface 
magnetization exponent reads 

^^ = 2 hr2 • 

The surface magnetization exponent is shown as a function of r/s in figure 6. The 

minimum value (3^ = 1/2, corresponding to the homogeneous system and to the bond 
problem behaviour as well, is reached when r = s. 



9. Conclusion 

We have presented a comparative study of the influence of bond and site aperiodicities 
on the critical behaviour at a second order phase transition. One dimensional 
sequences generated through substitution, corresponding to an uniaxial aperiodic 
modulation, have been discussed. The extension of the Harris criterion to the site 
problem requires the knowledge of the eigenvalues of a substitution matrix which 
is linked to the distribution of pairs of successive letters along the sequence. For 
sequences with IA2I > 1, the relevance-irrelevance criterion is the same in the bond 
and site problems. When IA2I < 1, the aperiodic perturbation may become more 
dangerous in the site problem, depending on the form of the substitution. 

Exact results have been obtained for the surface magnetization of layered Ising 
aperiodic systems. The period-doubling sequence leads to a marginal perturbation for 
both problems, but with different nonuniversal exponents. The Thue-Morse aperiodic 
modulation, which is irrelevant for the bond problem, becomes marginal for the site 
problem, where the surface magnetic exponent is nonuniversal. 

Among the possible extensions of this work, one may mention the treatment of 
substitutions with more than two letters and, more interesting, the study of higher- 
dimensional aperiodic perturbations in isotropic or anisotropic critical systems. 
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